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I.  Surtinary  and  Introduction 


linbedded  in  additive  Gaussian  noise,  an  unknown  acoustic  signal 

0(t)  is  recurring  over  and  over  again.  The  signal  consists  of  a 

specified  number  of  amplitude  modulated  pulses  of  width  T  sec's.  That 
n 

is,  0(t)  ■  Z  0  P[t-(i-l)T]  where  n  is  known  and  P(t)  - 

i-I  ^ 

We  do  not  know  these  amplitudes  0^,  but  suppose  Z  0^  and  the  average  rate 

of  recurrence  are  known.  Assume  that  the  recurrence  times  for  the  signal 

are  purely  random.  We  wish  to  determine  these  recurrence  times,  i.e.  , 

to  detect  when  the  signal  appears  in  the  noise. 

Let  X(t)  denote  the  signal  plus  noise  process.  We  sample  X(t)  by 

taking  non-overlapping  discrete-time  records,  each  being  wT  sec's  in 

duration.  To  be  more  exact,  each  record  consists  of  w  successive 

observations  on  X(t)  where  the  observations  are  taken  T  sec's  apart. 

The  optimal  method  for  determining  whether  or  not  0(t)  is  present 

in  a  given  record  is  given  by  the  likelihood-ratio  test  (Mood  [h], 

Wainstein  and  Zubakov  [^J]).  If  the  signal-to-noise  ratio  is  low,  then 

the  likelihood-ratio  test  amounts  to  observing  whether  the  sum  of  all 

the  observations  in  the  record  exceeds  in  absolute  value  a  fixed 

n 

threshold,  provided  that  Z  0 

i-1 

results  are  given  for  the  case  where  Z  0^  ■  0. 

The  detection  problem  is  closely  connected  with  the  estimation  of 
0(t)  and  its  autocorrelation  ^  (T)  ■^©(t+t')  0(t)dt.  Hinich  [2] 
discusses  this  connection  and  using  the  tools  of  large-sample  theory, 
shows  that  the  optimal  estimator  of  ^  is  a  linear  combination  of  the 
sample  autocorrelation  and  the  square  of  tiie  sample  average  of  the  record 
again  given  low  signal-to-noise  ratio.  However  the  estimation  of  0 


-  /  0(t)dt  is  non-zero.  However 


ri  if  0  <  t  <  T 
\0  otherwise. 


is  more  involved. 


Although  this  work  should  be  of  some  use  in  the  problem  of  active 
and  passive  sonar  detection  of  submarines,  the  initial  motivation  was 
from  an  investigation  into  an  underwater  communication  system  based 
upon  acoustic  pulse  position  modulation,  AFM.  To  Illustrate  this  system 
suppose  submarine  A  is  sending  a  message  to  submarine  B  by  APM.  Sub  A 
repeats  an  acoustic  signal  W(t)  as  defined  above.  The  intervals  between 
the  repetitions  of  y(t)  contain  the  information  which  A  is  sending. 

The  message  is  coded  so  that  these  intervals  seem  to  be  purely  random  in 
length.  However,  the  average  time  between  recurrences  of  the  signal  is 
at  least  an  order  of  magnitude  longer  than  the  duration  of  the  signal, 
nT  sec's. 

Sub  B  knows  n,T,  and  the  pseudo-random  code  for  the  intervals 
between  recurrences.  If  Sub  B  knows  the  0^,  it  could  detect  the  times 
of  occurrence  by  matched  filtering  (Wainstein  and  Zubakov  [5])» 
Unfortunately  the  medium  often  distorts  and  delays  the  pulses,  and  thus 
Sub  B  then  doesn't  know  the  shape  of  the  transmitted  signal. 

To  concluoe,  let  us  outline  this  paper.  In  Section  II  we  give  a 
formal  statement  of  the  problem  posed  above.  In  Section  111  we  present 
(as  a  Bayes  strategy)  the  likelihood-ratio  test  for  the  detection  of 
the  recurrence  of  the  waveform.  In  Section  IV  we  give  an  approximate 

likelihood-ratio  test  for  the  detection  of  the  waveform  when  Z  0.  “0. 

1 

Numerical  examples  are  given  for  both  the  case  where  Z  0^  /  0  and  when 
it  is  zero. 

I J .  Statement  of  Problem 

We  shall  develop  a  formal  statement  of  the  problem.  First  let  us 
discuss  it  informally. 


Wo  observe  n  process  X{t)  which  consistB  of  fi  rsndomly  occurring? 
unknown  signal  j.lus  noise.  The  noise  process  N(t)  is  assumed  to  be 
stationary  and  Gaum  an  with  mean  zero  and  known  covariance.  Without 
loss  of  tenoiality  wo  nay  normalize  so  tiiat  the  noise  has  variance 
°N(t)  ■  Otherwise  define  X*(t)  -X(t)/oj^^^^. 

The  waveform  0(t)  has  known  length  and  can  be  represented  as  a 
step  function  as  is  illustrated  in  Figure  1  for  n  ■  2. 


The  parameter  T  is  the  pulse  width  and  n  is  the  number  of  pulses. 

Thus  nT  is  the  duration  of  the  signal.  The  vector  of  pulse  amplitudes 
9'  *■  (W,  ....  ,9  )  is  unknown ,  but  assume  we  know  9,  ■»  ♦  ♦  •  •  9  . 

We  assume  that  the  time  intervals  between  repetitions  are  large 


compared  to  the  length  of  the  waveform.  Let  y  be  the  rate  of  repetitions 
of  9(t).  Then  ~^T  seconds  is  the  average  time  between  waveforms. 


ind  thus,  is  small. 


n 


1 

? 


We  assume  that  It  9  ll  “  (  Z  9.)  is  small  compared  to  the 

i“l  ^ 

variance  of  ti,e  noise,  we  can  state  this  in  terms  of  the  signal-to- 


df 


noise  ratio  of  X(t).  By  definition  R  •  ^  L  0./  /  Sij(f) 

W  ^  '^-00  N 

where  S  (f)  is  the  spectral  density  of  the  noise,  N(t).  But 
N 

2^j(f)  df  *  tlN^(t)  ■  ^  normalization  of  N.  Thus, 

-  -cx)  ‘  n  p  P  ^ 

R  =  Z  0.  •  ^  11  9  il  .  Since  ^  and  11  9  ll  are  small,  R  is  low. 

I'l 


Now  we  will  discuss  the  sampling  procedure.  For  fixed  integer  w, 

we  will  take  a  finite  group  of  w  successive  discrete  observations  on 

X(t);  the  successive  observations  being  T  second  apart.  That  is  for 

each  t,  such  that  t,  <  t.  <  *••  <  t  ,  we  observe 
i  12  m* 

X(t.  4  T),  X(t^  4  2T),  ...,  X(t^  4  wT). 

The  intervals  between  the  ^^i4l  “  ^i^  substantially 

greater  than  wT  and  nT  (the  duration  of  the  waveform).  Let 

X^^^  -  (X(t^  4  T),  ...,  X(t^  4  wT))  i-l,  ...,  m. 

We  thus  have  a  sample  of  ra  vector  observations  on  a  w-dimensional 
random  variable. 

Tills  sampling  scheme  may  be  regarded  as  opjening  a  sequence  of 
windows  of  width  wT  seconds,  through  which  we  observe  the  process  at 
ra  different  stages  of  time. 

There  are  several  different  possibilitj es  which  can  occur  when  a 
window  is  opened.  There  may  be  no  part  of  the  waveform  present  during 
the  wT  second  "look"  at  X(t).  In  that  case  we  observe  only  noise. 
However,  the  window  may  open  just  as  the  front  part  of  the  waveform  is 
"visible".  In  that  case  only  the  head  of  the  vector  ©  plus  noise  is 
observed.  Similarly,  we  might  observe  only  the  tail  of  0  plus  noise, 
or  perhaps  the  middle  of  0  plus  noise.  Incidentally  let  us  suppose 
that  wT  «  ^  T  since  we  wish  to  exclude  the  possibility  of  catching 
two  successive  waveforms  in  the  window. 

Since  the  window  has  w  components  and  0  has  n  components, 
there  are  n4w-l  ways  of  catching  part  of  0  along  with  the  noise. 

We  can  represent  these  n+w-l  possibilities  for  0  in  the  window  by 


defining 

(1) 


-5- 


) 


where  ©,  -  0  if  k  <  0  or  k  >  n^l.  For  example  (S  0) '  ■  (O  0,  ...,0) 

k  -  -  n-1  n 

and  (S  “  (0*  Of  •••t  There  is  no  way  of  knowing  in  advance 

which  case  is  occurring,  riach  of  the  n^w-1  possibilities  may  be  regarded 
as  equally  likely.  Figure  2  gives  an  exariple  for  w=3f  n*2  where  the 
noise  has  been  removed. 


Figure  2,  Example  of  Sampling  System  with  Noise  Kemoved  and  n“2,  w«3* 


'.Ve  will  catch  some  part  of  ©  i.i  the  i^^  window  if  and  only  if 

©(t)  commences  at  time  t  :  -nT  -ft,  ♦T<t  '^t.  ♦  wT,  Thus,  the  probabil- 

o  i  o  i 

ity  of  this  event  is  approximately  (n+w-1)  ^  where  is  the  recurrence 
rate.  The  probability  for  observing  a  specific  one  of  the  n^w-1  possibil¬ 
ities  is  simply  ^  , 

‘iince  the  distances  between  windows  are  greater  than  nT,  a  single 

waveform  cannot  appear  in  two  successive  windows.  Moreover,  suppose  that 

for  some  »  bN(t*  T)N(t)  =  0  forX  >‘T  ,  Then  if  we  take  the  windows 

o  o 

further  apart  than  'X' ^  seconds,  the  X^^^’s  are  independent.  The  above 
restraint  on  the  covariance  of  the  noise  holds  approximately  for  many 
colored  noise  processes  which  occur  in  applications.  Of  course,  for  white 


-6- 


noise,  'T'  -  0. 

We  then  can  sum  up  this  discussion  with  a  formal  statistical  state¬ 
ment  of  an  idealized  version  of  the  problem: 

W'e  have  m  independent  w-dimensional  random  vectors 
\  each  identically  distributed  as  X  where 

N  ♦  S  ^^2^0  with  probability 

N  ♦  S  H  n  .. 

-w+2 


(2) 


n-2 

N  >  S  ,© 
n-1 


l-(n^w-l)  / 


\  N  «  M 

and  SO  is  defined  in  (1)* 

J 

The  vector  random  variable  N  has  a  w-dimensional  multivariate  normal 
distribution  with  mean  zero  and  known,  non-singular  covariance  matrix  E. 

W'e  express  this  by  (N)  =  7/(0,  E)  where  ^  {n}  is  the  distribution  function 
of  the  random  variable  N, 


Suppose  that  E  is  either  zero  or  it  is  a  known  constant  of  an  order 

of  magnitude  greater  than  ||o||  ^  ~  S  0^.  Assuming  that  ^  and  ||  0  ||  are 

1  ( i ) 

.small  we  desire  to  test  for  each  record  X'  which  of  the  two  following 


hypotheses  are  true: 

(3) 


H  : 
o 


K^:  X^^^  » 


S^O  for  some  J  =  -w+1,,,,,  n-1. 


In  other  words,  we  wish  to  determine  which  records  contain  part  of  the  wave- 


foimi  along  with  noise 
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III,  Statistical  Detection  of  Wavefom  for  0^  ^  /  0 

In  tnis  section  we  will  derive  the  optljnal  (Bayeo)  test  of  the  two 

hypothesis  formulated  in  Section  II, 

Suppose  we  have  a  random  variable  X  from  a  population  with  density 

function  f(x  |  (|l).  Suppose  we  wish  to  test  against  (|  ■  (jl^^ 

where  (j)^  and  (|>^  are  completely  specified.  Let  h^  and  h^  denote  the 

a-priori  probabilities  for  (|)^  and  respectively,  i.e, 

p(d)  true}  =  h 
^~o  ^  o 

U) 

p((j»^  true)  “  -  1-h^ 

Let  L^  denote  the  loss  incurred  in  rejecting  when  in  fact  (|)  ■ 

i.e,  the  false  alarm  cost.  Let  L^,  be  the  loss  in  accepting  when 

i.e,  the  cost  in  missing  the  signal.  The  Bayes  strategy  is  a 

function  of  the  observed  random  variable  which  chooses  either  H  or  H-  in 

o  1 

such  a  way  as  to  minimize  the  expected  loss.  In  Chapter  12,  Mood  [4] 
shows  that  the  Bayes  strategy  is  to  reject  (accept  (|)  ■  (|)^)  if 


(5) 


,,  , 


>  k 


h  L 
o  o 


The  random  variable  X  and  the  parameter  can  be  vectors.  We  call 
\(x  )  the  liklihood-ratio. 

We  will  now  deal  with  the  densities  of  interest  in  this  work.  However, 
to  facilitate  the  algebra  we  make  the  one-to-one  transformation, 

(6)  z  =  X 

where  is  the  inverse  of  Z,  the  covariance  matrix  of  the  Gaussian  noise 


vector  N 
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From  (2)  we  have 


N*  *  Z~  SO  with  probability  for  each 

«J 

J  -  -w+1, . . .,n-l 


^  N*  with  probability  ^l-(n^w-  ". ) 
where  ^{N*}  ■  Tv  (O,  Z  ^). 

If  the  noise  N(t)  is  white,  then  EN(t^)  ”  0  for  t^  /  t^. 


Thus  Z  is  the  identity  matrix  since  we  made  o! 


N(t) 


1,  Thus  Z  -  X 


To  restate  the  hypotheses  and  in  terms  of  Z,  from  (3)  v^e  have 

H  :  Z  ”  N*  with  a-priori  h  -  l-(n^w-l)^ 

°  -1 

Z  ■  N'^  ♦  Z  SjO  for  some  J 

Thus  says  that  only  noise  is  present  in  the  window  and  says  that 
some  translation  of  ©  is  present. 

Let  f (z 1  0)  be  the  density  of  Z  given  parameterized  by  the 
waveform  vector  0.  Therefore  from  (7)  we  have 
(9)  f(z|Q)  =  — Z  n(7,|z"^5  0,  Z"^) 

J=-w*l 


where 


•  ?  ,  -  7  '7 

n(z|(p,i;)  =  (ZTT)  ^  |c|  « 


is  a  w-dimensional  normal  density  with  mean  (p  and  covariance  matrix  C, 

The  density  of  Z  given  is  simply 
f (z  I  0)  =  n(z  I  0,Z"^) 

Notice  that  f(z  |  0)  is  a  convex  combination  of  multivariate  normal 
densities,  but  it  is  not  in  general  multivariate  norrmil  itself.  We  shall 
handle  it  by  muaking  Taylor  series  approximations  with  0  in  the  neighbor¬ 


hood  of  zero. 


-9- 


Dy  expanding,  in  Taylor  series  about  0  =  0  we  have 


(10) 


n(z  I  0,  l"  ) 


♦  iia)ii^  K*(z,(p  ) 


where  the  c's  are  constants  and 


Z-i 

♦.  .  .♦  t* 

z 

K’^(z,  (fj  ) 

,  ± 
d  e 

1 

w 

w 

for  some  t*  >  0  ard  d  >  0.  Thus  b^[K’*(Z,(|)  )]^  exists  and  is  bounded 
by  some  namber  independent  of  (p  for  each  r  >  0.  With  the  notation 


(11) 


p  /  .  .  ,  ^  ij  ik  Jk 

G^(z  i,j,k)  -  z.z^zj^-o  z^-o'"  z^ 


E^G^(z I i, J,k)  =  0  for  all  i,J,k 


From  par,e  39  of  Anderson  [1], 
o  I  j  ‘k^l 


, . ,  iJ  k  /  ^  ik  j/  i  /  Jk 

h  =  o  '’o  *  a  o  ♦  o  o'' 


K  Z.Z,  =  o^*^. 
o  1  J 

Moreover  all  odd  moments  of  the  are  zero.  We  then  have  the  following 


orthogonality  reiationchips ; 


ik. 


K  Z.  (Z.Z,  -0  )  =  0 

o  1  J  k 

K^Z.G^(Z  J,k,  0-0 
r;^(Z.Zj-o"'^)G^(Z  1  k,A  ,m)  = 


for  all  i,j,k 
for  all  i,  J,k, 

0  for  all  i,  J,k,^'  ,m. 
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}’uttinf;(p=  OjO  in  (10)  and  summinf;,  we  have  from  (5)  and  (8)  that 
the  likJihood  ratio  K(z  j 0)  =  ~f |u j 


(12) 


f(z 

0} 

f(z 

0) 

J  cl 


♦  Z  G  (z|i,j,k)0(||o|P) 

i.  J,k 


>  K(z,«)  0(11011**)]. 

E^[K(z,0)J  exists  and  is  bounded  by  sone  number  independent  of  0 
for  each  r  >  0.  Moreover  the  0(  jjojp)  and  0(  |j  0  |p)  terms  are 
functions  of  0  v;hich  do  not  involve  z. 

Now  define  the  discrete  autocorrelation  vector  V  *  =  (”7^^, . . . , 
and  the  time-averar.e  i/  by: 

i  («?  *  4  4) '  i  il« 


(13) 


1  2  '1  2  n'  2 

>'  =Q^  *0  Q 

2  12  23  n-1  n 

V  n  =  0,0-,  ♦  0.,0,  ♦...♦«  ^0 

T  3  13  2  U  n-2  n 


T  n  i  r 


n 


(14)  V  Z  0.  (DC  value  of  0). 

°  i  =  l  ^ 

Applying’,  (1),  (13),  and  (14)  we  have 


(15) 


n-l 

;  D  u  =l0  1 
. .1  'O  w 


J  =  -v^l 

w 


where  1'  l,...,Jj,  nnci  for  any  syiijaetric  w  x  w  matrix  A  ■  (a.  .) 

w  1 J 
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(16) 


^  (.s  g)‘  A(L-  U)  =  (  L  a 

'  J=-wl  J  i=l  ^  i-1  ^ 


•  •  • 


w-n+1 


i-1 


where  is  understood  to  be  7,ero  when  i  or  J  is.  greater  than  w  or 

less  than  1.  Tlius,  if  w  <  n,  the  coefficients  of 
vanish  and 


w 

^  i;  (b  g)'  A(b  e)  -  (  t.  a  jV 
^  j  j  j  i=i 


Applying  (13)  and  (16)  to  (IP)  we  have  the  following  result, 
Le.T.’na  1,  Define  Y(z)'  =  (Y^  (z ) , , , .  ,Y^(z ) )  by 


(17) 


Y,(.) 


w-k  +  1 
L 

i  =  l 


i,i+k-lv 

j  k-1, , , ,  ,n. 


Then  the  likl ihood-ratio 

(18)  |»)  =  =  1  .  ia':„)V„  •  r(z)'r 

-  I  Go(z(i,J,k)0(||y||^)  ♦  k(z,«)0(  ||o  K^) 

ij,k  ^ 

It  IS  understood  that  z  .  ■'=  0  if  i  o  or  i  >  w»l,  and  thus  if  w  <  n 

1  -  - 


Y  , (z)  =  ...  =  Y  (z)  =  0. 
w+1  n 

(-uppose  we  know  th'“  waveform  vector  y.  Then  from  (5)  and  (8)  in 
;  ection  11,  the  Hayes  test  between  (only  noise)  and  (some  part  of 
the  waveform  pres('nt)  is: 

Reject  if  \(z|  y)  >  k. 

flowever  supro.*  •>  we  know  onlyl/)^  =  10^  but  we  use  an  alternative  test 
whicfi  says  rs'  'ect  H  if 


-L?- 


(19a)  K^Czle)  =  1  ♦  T  (z'l  )V  >  k. 

'  '  n^vr-1  w' 

approximation  of  From  Hinich  [2],  we  have 

--1  .  ^  '^o 


\*  is  the  linear 


E_  Z'l 
Q  w 


(!'  Z 
w 


n-'-w-l 


Var_  Z'l  «=  1'  Z~h  *  -  0(||e||^). 

6  w  w  w  n^w-1  M  II 

Thus  we  can  rewrite  (19a)  to  say  re.lect  if 
(19b)  (!'  Z"^  1  )  Z'  1  >  (k-1). 

W  W  W  '\jj  ^  ' 

Under  H  ,  Y  =  (!'  Z~^  1  )  Z'l  is  a  normal  random  variable 

o  w  w  w 

with  mean  zero  and  variance  one.  Under  Y  is  again  normal  with 
variance  one  but  with  mean 

(20)  (1'  Z^^  1  )  . 

'  w  w'  n+w-1 

Let  a  and  g*  be  the  probabilities  for  rejecting  when  in  fact 

it  is  true  (false-alam)  for  the  \  and  X*  -  tests  respectively.  Let 

3  and  3"^  be  the  probabilities  of  accepting  h^  when  is  true,  i.e, 

the  probabilities  of  missing  the  waveform.  Then  the  expected  loss 

£(L)  using  \  is, 

(21a)  E(L)  =  gL  ♦  3L, 

o  X 

where  L  is  the  false-alanr.  cost  and  L,  is  the  cost  of  missing  0. 
o  1  “ 

Similarly  the  expected  loss  E(L*)  using  \*  is. 


(21b)  £(!■")  =  g’L^  ♦  3*1^. 

Since  the  \  -  test  is  Bayes,  E(L)  <  E(L*),  But  in  the  Appendix  we  prove 
Theorem  I:  The  -  test  (19b)  between  and  has  expected  loss 

(22)  S(L^)  <  E(L)  ML  ♦  L  )  - 2 - 5-  L  V'D  0(  ||  0  j(^)  ]  ♦  0(  ^  ) 

°  (w'^n-l)  " 
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whero  D  is  the  n  *  n  matrix 
(23)  D  -  Y(Z)'. 

!' Example  1:  Suppose  the  noise  is  white.  Thus 

};  =  1  nnd  X  =  Z 


w 


'Hien  from  (17)  we  have  the  lutocorrelations 
(X)  =  (X‘^  -  Ij  (x||  -  1) 

Y,,  (X)  =  X-X.,  ♦...♦X  -  X 
2  ^  12  w-1  w 


t 


{2U) 


Yn(X) 


XX 
1  n 


X 


w-n^l 


X 

w 


Luprose 


Thus  from  (13)  and  (1/*) 


(25) 


2  ^2*  n— 1) 

n  ■ 


'V'  o  =  - 

het  w  =  m.  Thus  fron  (23)  and  (24) 


Thus  from  (25)  and  (26) 

y  ■  L-v  s 


Now  let  L  =  1  and  L,  =  3.  This  means  that  a  miss  is  three  times 
o  1 

as  costly  as  a  false-alarm.  This  is  realistic  if  there  is  a-priorl 


knowle(i/’e  of  the  pattern  of  repetitions  of  0,  allowing  us  to  reduce  the 


-u- 


false-alarm  rate. 

Moreover,  assume  that  l/a  *  5n,  i.e. 
occurences  of  6  is  five  times  as  long  as  6 
alone  in  the  window  is  Thus  we  have 


the  average  distance  between 
.  Then  the  probability  of  noise 


1  L  , 

1-na  o  _  4 

na  3 


Now  let  n  «  w  =  9.  From  (19b)  we  reject  if 


(27)  ^  L 

1  I  o 

1  ^  I 

Under  H  ,  Y  »  —  f  X  has  density  n(y  |o,  1  ), 

O  1  ^ 

(28a)  a*  -  P^{Y  >  1.89|h^}  -  .029 

Under  Y  has  density  n(y|d,  1)  where 

d  “  — ■'III  =  .18  from  (20).  Thus 
n^w-1  yo  ' 

(28b)  0*  -  {Y  <  1.89|h^}  -  .956 


From  (21b) 

E(L*)  -  2.89 

From  (22)  and  since 7, 

U 

(29)  2.65  <  E(L)  <  2.89 

From  (28b)  we  see  that  the  probability  of  a  miss  is  very  high. 

This  is  because  the  detectability — d  is  small.  But  this  is  the  best  we 
can  do  given  this  sampling  scheme.  The  signal-to-noise  ratio  is  small. 
If  we  used  matched-filtering  with  0  as  the  filter  function,  we  could 
not  do  better, 

1 

The  test  statistic  E  is  easy  to  compute.  If  we  know 
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yj  we  can  easily  carry  out  the  X’*'  -  test.  For  small  signal-to-noise 
ratio  it  is  optimal.  It  should  still  be  reasonably  good  even  if  the 
signal-to-noise  ratio  is  not  too  small,  especially  if  we  include  the  second 
order  terms,  which  involve  the  waveform  autocorrelation  terms,  ■'IjJ  and 
the  autocorrelations  y^(Z). 
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IV,  Detection  of  VJavefom  vdthV  “  0. 

*  o 

If  TjJ  o  =  0»  the  \*  -  test  is  no  good.  But  suppose  we  know  the  discrete 

autocorrelation  vector  .  Furthermore,  assume  the  noise  is  white.  Thus 

E  »  I  and  Z  =  X,  To  test  between  H  and  H,,  let  us  use  the  alternative  test 

o  1’ 

which  says  re.ject  if 

(30a)  (zIq)  °  1  >  — Y(X)V>k, 

where  the  vector  Y(X)  is  given  by  (24)  and  E^YY'  =  D.  X.**  is  the  second- 
order  aporoximation  of  the  likelihood-ratio  X. 

From  Cochran's  theorem  and  the  central  limit  theorem  (Lindgren  [3]  )» 
und*»r  H 

o 

(31)  Q  “  (-V  "  Y  (X)y 

J-1  ^  ^ 

has  a  normal  distribution  with  mean  zero  and  variance  one  as  v-*CP  . 

Thus  we  can  rewrite  (30a)  to  say  reject  if 

(30b)  Q>  (k-1)  (n-'-w-l)  (^  ' 

Under  Q  is  again  approximately  normal  with  variance  one,  for  large  w, 
but  with  mean 

n>w-l  (  V 

since  frxam  Hinich  [2], 

’'e  0(||e||2). 

Vie  thus  can  prove  in  a  manner  similar  to  that  used  for  Theorem  1, 

Theorem  2;  For  small  ©  ,  the  X""*  -  test  (30b) 

has  the  property  that 

E(L)  =  Ed’’")  = 
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where  u*'*  Is  the  fnlse-nlarm  probaDility  for  the  \**  -  test  and  3**  is 
the  probability  of  a  miss.  In  this  case  we  do  not  specify  the  error  terms. 
However,  by  analo/^  to  Theorem  1,  the  \**  —  test  should  be  nearly  optimal 
for  reasonable  larp.e  n, 

Kxample  2:  For  w  >  n,  from  (23)  and  (24) 


Let  w  =  8  and  n  =  4  with 


«'  =  (1,  -1,  -1,  1) 


From  (13)  and  (14) 


.  0  andV  -  i  (2,  -1,  -2,  1). 


Thus  from  (33) 

- ,  '  n, ,  .  loo 

y  DV  =  -yg— . 

Let  k  -  1  =  y  as  in  Lxample  1,  From  (30b)  and  (31)  we  reject  if 


Q  >  1,47  where, 


Q  =  i  [2  £  -  1)  -  r  -  2  I  Xj  X..J  .  r  x^  x^.j  ]. 


Under  H  ,  Q  has  density  n(q  0,  1), 


(35a)  =  P  {Q  >  1.47  I  H  }  ■»  .071 

r  '  o 

Under  H^,  w  has  density  n(q|d,  1)  where 

d  =  - ^ —  ( V  D'V')  =  .23  from  (32).  Thus 

n+w-1  '  ' 

(35h)  b"’  =  ^  1.47  H^)  =  .893 

Thus  fror.;  Tiieorem  2, 

:'.(L’*)  =  2.75 
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To  conclude,  we  see  from  (31)  and  again  in  the  example  with  (34)  that 

the  test  statistic  is  simply  a  cross  correlation  between  the  sample  discrete 

autocorrelation  vector  Y(X)  and  the  discrete  autocorrelation  vector  .  This 

cross  correlation  is  fairly  easy  to  compute  provided  we  know  the'V/^.  Hinich 

[2]  gives  the  optimal  estimators  of  the  V  based  upon  Y(Z)  and  Z'l  . 

1  w 

In  addition  we  assumed  white  noise  for  this  section.  The  case  of  a 
general  covariance  matrix  I  is  not  a  difficult  extension,  it  just  complicates 
the  algebra. 


Appendix 


(Al) 


Proof  of  Theorem  1;  From  (10) 

a  -  P  {\  >  klH  }  >  P  (a**  >  k  ♦  -,lx  -  \*|<  -|h  } 
r  '  o''  -  r  '  nl  o 

-  P  {\*  >  k  ♦  -|h  }  -  P  (\*  >  k  ♦  -,K  -  \*|>  -|h  } 

r  ^  n'  O'*  r  ^  n*l  I  n|  o'* 

>  P  {(!'  1  )  L  Z.  >  k‘|H  }  _  p  {|\  -  \*|>  -^Ih  } 

-  r  ^  w  i  -  I  O''  rM  I-  nl  o'* 

where, 


k'  - 


(1 


n^w-1  I  -1  /.  ,  .  1\ 

775777^  "t'  o 
w  w 


Now  a*  -  P  {\*  >  klH  } 
r  I  o' 

where  k''  ■  _ n»w-l _ ^  (k-1), 

1/2  ^  ® 

(!•  r  ^  1 .) 


Since  (1 


w  w^ 

,  £-1  ,  -1/2  “ 


w 


Z,  under  H  is  normal  with  mean  zero  and 
i  o 


(A?) 


variance  one,  and  n  is  large, 

{(1;  1^)  r  >  k-}  -  a*  *  o  (  i  ). 

Moreover  by  the  Tchebysheff  inequality.  Mood  [/*), 

2 


f'r{  I  ^  -  ‘*1^  k  l«o>  5  ,  "  o  t t'DY  *  0  (  )|  e  II  “)] 

(w+n-1) 


since  (\  -  \*)  =  0  and  from  (18)_^ 

—  [  v  Dy/  >  o(||  e||^) 


(A3) 


E  (\  -  K*)'^  -  - ^ 

^  (w+n-1) 


Thus  from  (Al)  and  (A2), 

N  *  2 

a  >  a’  -  n 


[yj'  0  (lie  11^)  ]  ♦  o  (  i  ) 


(w>n-l) 

Similarly  for  3  =  P^  {\  <  k|H^} 


2  L 

(A/.)  (!>()•-  - - 5—  tv '  uy  ‘  0  ( ||e  II  ) 

(w^n-1) 

Since  from  (21b),  E(L*)  ■ 
we  have  the  desired  result. 
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